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ABSTRACT 


The goal of our research is to investigate the collective modes of the anyons localized in 2D 
parabolic well and get the right exact expression for their collective mode frequencies. A study of the 
collective motion of atomic gases, localized in the harmonic trap, belongs to class of actual and 
interesting problems of physics of ultra-cold atomic and molecular systems. A topology of 2D 
systems allows to exist the particles, whoes statistics may be orbitrary between bosons and fermions, 
therefore they call anyons. And these anyons are described with parameter and may be determined 
in the interval between ones for bosons and fermions. The one of intriguing problem of ultra-cold 
atomic gases is a study of the role of the anyon statistics to the system centre mass mode frequencies. 


KEYWORDS: 2D, Anisotropic. 


1. INTRODUCTION 


It is well-known that all elementary particles fall into one of two possible categories - bosons and 
fermions, depending on whether they obey the Bose-Einstein or the Fermi-Dirac statistics 
respectively. These particle are at least in 3-dimensional space-time. However in two space 
dimensions we do not have only bosons and fermions, but also particles with any statistics in 
between. These particles are called anyons and are the subject of this work. 


Certainly, it is unusual feature of anyons that they arise only in two-dimensional systems and it is 
hard to imagine for both physicists working at totally different field of the physics and people far 
away from science these amazing particles. However, these particles are not simply topological 
fantasies or objects of purely mathematical interest; on the contrary they might play an important role 
in certain physical phenomena of the real world. Of course, since we are living in at least three space 
dimensions where particles can be only bosons or fermions, anyons are not real particles. However 
there exist certain condensed-matter systems (for example thin layers at the interface between 
different semiconductors) that can be regarded effectively as two-dimensional. Their localized 
excitations (if they exist) are quasi-particles subject to the rules of a two-dimensional world. It is 
these quasiparticles that may be anyons and may be observed in certain cases. For example the 
collective excitations above the ground state of systems exhibiting the fractional quantum Hall effect 
(for a review see (Prange and Girvin 1990)) have been identified as localized quasi-particles of 
fractional charge (Laughlin 1983), fractional spin and fractional statistics (Arovas et al. 1984; 
Halperin 1984), and thus they can be regarded as anyons. Furthermore, anyons are conjectured to 
play a role also in the theory of high temperature superconductivity (Chen et al. 1989), even though 
in this case no conclusive word can be said at the moment (Lyons et al. 1990; Kiefl et al. 1990; 
Spielman et al. 1990). 


Such anyonic particles are becoming of increasing importance in condensed matter physics and 
quantum computation. They may play an essential role for describing the fractional quantum Hall 
effect, high-temperature superconductivity, and the physics of topological insulators and 
superconductors. Moreover, anyons as unusual quasiparticles with properties of its statistics are 
adequate tool for implementing a topological quantum computer. All of the mentioned categories in 
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physics are the active issues of the physicists throughout the world. Certainly, all of these emphasize 
the importance of investigating the quantum anyon systems and their collective motion in the 
harmonic trap that has been done in this scientific work. 


Most of the great interest that anyons have attracted in the past few years derives from the 
(unexpected) applications of these ideas to certain two-dimesional condensed matter systems, most 
notably those exhibiting the fractional quantum Hall effect (see for instance (Prange and Girvin 
1990)). In this case a series of new states of matter emerge as incompressible quantum liquids 
(Laughlin 1983) around which the low-energy excitations are localized quasi-particles with unusual 
fractional quantum numbers, i.e. anyons. Furthermore, it is also very likely that anyonic excitations 
with fractional statistics exist in films of liquid 3He in the A-phase (Volovik and Yakovenko 1989). 
The application of anyons to the theory of high temperature superconductivity has also been 
considered quite extensively (for reviews see (Wilczek 1990; Lykken et al. 1991)), but their actual 
relevance in this context is quite controversial and doubtful. 


Since experimentally first ultracold atoms have been realized in harmonic potentials, the goal of our 
work will be the consideration of the collective motion of anyons in the 2D harmonic trap. 


The present paper is organized as follows. We start with introducing Hamiltonian of the anyons in 
2D parabolic harmonic well in section2. Then, cumulant method is introduced in section3. In 
sections 4 and 5 some calculations have been given by utilizing this method. Next section is devoted 
for deriving the harmonic oscillator equation for the centre-of-mass - the main result of our work. 
Finally, at the end the conclusion is presented. 


2. HAMILTONIAN OF ANYONS TRAPPED IN 2D ANISOTROPIC HARMONIC 
POTENTIAL 


In this section we describe the Hamiltonian of anyons, localized in the 2D anisotropic trap, which 
expression will be taken from the paper [26] and, following to the paper of Ghost and Sinha [27], we 
write this system Langrangian. 


The Hamiltonian of the gas of N anyons with mass m and charge e, confined in 2D parabolic well, 
iS: 


N N 

a) 1 — —>- 2 m 

H= —) (Dx + Ay(7%)) + > =F (wei + wiyp). 
k=1 k=1 


Here 7, and px represent the position and momentum operators of the k th anyon in 2D space 
dimension, 


Cz X TR; 


Ay Ge) = hv) 
j#k Irj | 

is the anyon gauge vector potential [28], 7; =7 —7 and é; is the unit vector normal to the 2D 
plane. In the expression for vector potential A,(7%),v is the anyon factor and hereafter we assume 
that 0 < v < 1, which means the variation of the anyon factor between bosonic and fermionic limits 
of anyons. 


Our interest is the solution of the Schrédinger equation 


Let us consider first the term in the Hamiltonian A, containing only the anyon vector potential 
A, (rg). In the bosonic representation of anyons we take the system wave function in the form 
[29,30] : 
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Here and above R = {7{,T>, -.»M} is the configuration space of the N anyons. The product in the 
right hand side of this equation is the Jastrow-type wave function. It describes the short distance 
correlations between two particles due to anyonic (fermionic) statistics interaction. 


By substituting the wave function of this form into Schrédinger equation (3) without the harmonic 
potential term, we obtain the equation: 


aW,.(R, t th a z 
i a) = (A, + A,)¥7(R,t) 
where 
N go 28 
A ys h7 Ay h2v Vkj : Ve 
si - Nee 4 (2 
k=1 ae a j#k [7 
and 
N 4 
nA h > 5 — A r " Tr; 7 
fy =-i=) Ay Ge) We +) as 
=r jek les 
ing 


As in the paper [27] of Ghosh and Sinha, by introduc- 


h 
Ao = 
MW 


where Wo = ,/W Wy, we make dimensionless the length quantities and denote them by tilde sign. 


We express the energy quantities in the Hamiltonian (1) in the units of Aw 9. Then, for instance, the 
harmonic potential term will have the form: 


N N 

m hwo Pe aes ig 

zd, (uted + 08) =D (rai? + 30%") 
where A = W,;/Wy, Xz = Xx/Q9,Vx~ = Vx/Ao and parameter A is the anisotropic parameter for the 
harmonic potential. 


Now we make dimensionless Hamiltonians H, and H, 


x N h7A, h2v Tej : Vie 
A, = > = 


k=1 em si j#tk lel” 
N A x v 
=hiy ) +y) a ~ 
k=1 j#k 7.)| 


since A, = 07/ 0x2 + 07/ Oyz = A,/az. 
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Similarly 
N Rs 
AG 
Hy = —ihwo ». Ay (7) Vi tv ve) hs 
k=1 j#k [Fi | 
where 


with t = wot. 


At the end of this section, we emphasize that the wave function WP; (R, t) contains the configurational 


space of N anyons vector R. Therefore, it corresponds to many particle wave function of system. 
Previously, at the calculation of time variation of BEC, the wave function was a function of only one 
coordinate of condensate (see, for example, the paper [27]) and the solution of problem of BEC 
collective motions in the harmonic trap was essentially easier. 


3. CUMULANTS EQUATION OF MOTION METHOD 


For the description of above mentioned monope and quadrupole modes and also the oscillation of the 
centre of mass motion (the Kohn theorem), we use the cumulants equation of motion method [31, 
34]. According to this method, for the small amplitude oscillations, it is convenient to take the trial 


many body wave function ¥; (R, t) in the Gaussian form (we use notations, taken from Ref. [31], for 


variational parameters): 
N 
[ | e|- (sats) x 
exp|—(=> + tA, 
nos 244 


1 
(Xe — Xo)? + ix, 6, — (a + ia) (Ve — Yo)? + 
2 


iypC] 
Here, x, and y, are the x and y coordinate components of k-th particle, all variational parameters 
41,92, Az, Az and C,, Cz, and centre of mass components xg and yg are the time t dependent. 


> 1 
i ie aa) 


In order to derive the comulants equation of motion, we average over the Schrédinger equation (11) 
the weight fy y : 
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N N 

; i aya te i 

if | | dxedyeti, wit = ue dred Vefdy 
k=1 k=1 


N 
x Ps 1 
k=1 


Since the wave function VY, Eq. (12), is a Gaussian, no zero these averages are only for two central 
moments. Averages with fb = x, — Xo and i = Vx — Yo provide equations to find the centre of 
mass motion. And averages with ff = (x, —%0)*-—qz and f~ = (Ox —yo)? — 93 provide 
equations to find the widths motion. 


4. AVERAGE QUANTITIES FOR f, OF IDEAL GAS OF PARTICLES IN 2D 
ANISOTROPIC HARMONIC POTENTIAL 


For the ideal gas of particles in 2D anisotropic harmonic potential, we have an averaged Schrédinger 


equation: 
N 
if | | axa Se = | [] dx, dy efily * 


k=1 
N 

* 1 A we 2 1 2 
k=1 


Using the way, Ref. [31], of calculation of this equation integrals, we obtain: 


OW, 
| [| AX, AY, (X~ — Xo) Pr sae 


1 
v(t + iA, q?x% +iz= 9 ‘<,) 


and 
N 
Or 
| [| Ax, AY (VR - Yo) Pr ar 
Yo q2 
w(e + iA,g3¥y +i 9 £6) 
Then 
N N 
[ [] rece — 20)” Aer = 
k=1 k=1 
—N(iC, — 2C,A1q1) 
N N 
i | | Ax,A YK (Ve - Yo) Pr ». A,r = 
k=1 k=1 
—N(iC, — 2C,A2q3), 
and 
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| | AX AY, (X~ — Xo) Pr > Xr = 
Nxoqt 


N N 
| | | AXA (Ve —y)¥3). Her = 
k=1 k=1 
Ny043 
5. AVERAGE QUANTITIES FOR f,,, WITH ANYON PART OF HAMILTONIAN Ah. 
We start with the expression for the square of modulo of wave function Yr (R, t), Eq. (12), (for the 


simplicity, everywhere below, we omit signs tilte). It equals to 
N 
os aD 1 Caen 
wot = (2a) [Tool 
Nl = aaa, 11 qi 


_Yko 
q3 


where Xxpq = Xp — Xp aNd Veg = Ve — Yo- 
First, we need to calculate the integral in the average quantities for lem = Yxo, related to term 
Tei Vi 
Yeo t to: EE = 
cs 


=> 12 
Vio|Pr(R)| Zz 


2 2 1 
Xj t+ Vkj 


where 
eer 
Zy = —2X4jXxo oo + ia) + iC, X_j - 
245 


Aes 
2YKjVKO (a a iy) a 1C2VKj 
2 


The expression for this integral is: 


> N 
n=(2-) y > 
S 14192 ! 


{| | dAx,dy,ax;dy; X 
2 


k=1 j#k 
2 2 2 
Yro Xko YVko jo  Yjo 

a yn Z1EXP | a 

Xie + Vij q4 42 4 42 
We introduce new variables x,j = x, — x; and Yu; = Ye — y; then dx; = —dx,; and dy; = —dyx; 
and taking into account that 

2 


exp |-—> 
| qi 4 
Vij — 2 ej Veo + “4 
2 
q2 


2 
M0... a 
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the exponential function in the Eq. (24) will have the form: 


2 2 ois ; 
fe - 2Xko 20 Xkj 2X ej Xko 2 


aq 44 q 
Vij — nein 
q3 


Next, substituting expressions for Z,, Eq. (23), and exponential function, Eq. (26), in Eq. (24) for 
integral | a then, using the formula, Eq. (44), at the integration of / is over Ax, dy, dXzj AY, ;, we find 
that only the last term iCzy,; of Z, gives a non zero contribution into / ae And its expression is: 

Ty = iN(N — 1)C3. 


At the derivation of this expression for I? _» we have used the formulas: 


2 


+00 p+oo 
Xkj 41,2 B12 
dx, V_; =—— e 2 “kj 2 Vkj — 
J J y2 + y2 
—0o0 Y—0oo kj Vkj 


21 By 1 
J Py a, (a, /2 — B,/2) 
for a, > PB, and 
+00 +00 x2, 7 B 
kj —six? Shy? 
AX, ; AVR; >——e 20K 2K = 
fe I. xg + Vig 
27 ay 1 
Vf B, (Bi /2 — a/2) 
for B; > a, which can be obtained, using expressions Eqs. (30) - (32): 
+00 +00 x2, a B 
kj 21,2 Ely? 
AX, AVR; me 20K 2 Pk 
he P xij + Vig 
At the calculation of that integral over dx, j» we use the formula 3.466.2 of the book [32]: 
te x _ 242 us 1p 2 pR2 
l Un ese Be ree See [1 — erf(6y)] 
with [Re 6 > 0, argu < 2/4|] and erf(x) is the error function. 


At the calculation of obtained integral over the dy,;, we use the formula 6.289.2 of the same book 
[32] of LS. Gradshteyn and I.M. Ryzhik: 


+00 B 
erf(Bx)e(F?-#?)*" xdx = ——__, 
J we 2u(u? — B*) 
at [Re uw? > Re B?, |argu < 1/4|] 
In the analogous way, one may calculate the expression for integral 17, for the average quantities of 
fey = Xxo- It equals to expression: 


IZ = iN(N- 1). 
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6. AVERAGE QUANTITIES FOR Fey WITH HAMILTONIAN ee 
We calculate the integral in the average quantities for Ey = Yxo, related to term 
ios HS [Pr |? 
Veo rAy (7%) * VePr = VVK0 22, 
Xj + Vk j 


where 
i 
Zo = 2YKjXko 2q2 + tA, a ICV xj — 


1 
2XKjVKO (oe Ss iy) + 1C2X,; 


Again, the expression for this integral is: 


N 
1 2 
an ) dD Sip J axe Prilvex 
Z2 4192 Li Ls kaVK ‘J Vj 


Pe 2 
x2 . 
a Z2exp |-" . _ Yio ae 7 — 2A 
Xiej a Viej Gq &@ & 


We follow the procedure of calculation, discribed from Eq. (24) up to Eq. (29), except of substituting 
expressions for Zz, Eq. (35), and exponential function, Eq. (26), in Eq. (36) for integral eS then, 


using the formula, Eq. (44), at the integration of LF over Ax, dy,dX,j;AY_;, we find that only the 
term —iCyy,; of Zz gives a non zero contribution into [ om We obtain 

Iy, = —ivN(N — 1)C,. 
Analogously, we calculate the expression for integral I7, for the average quantities of les = Xxo. It 
equals to expression: 

Iz, = ivN(N — 1)C. 
We calculate the integral in the average quantities for fev = Yxo, related to the last term in the 
Hamiltonian 5. It is: 


Ay (Teej 
Sey BO, yt = 
k=1 j#k [re 
5 —VriXkj + XK kj 


ROC , 
2 2 2 2, Dp 
k=1 j#k l#k (xi; + yiij) ie + Vi) 


where 
2 
xX . 
Lexp = exp| “ie — Hie — Sp _ 28 
qt q3 qt q2 
qi 4 


Introducing variables x,; and y,;, we expressed the first exponential function in Eq. (40) in the form 
of Eq. (26). Now, we introduce variables x;, and y;; then 
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qi qq Ga 


will have a form: 


2 2 2 
as - 3Xico - 3Ykco _ Xk — 2Xk1Xk0 _ 


aq 4 q 
Yin — a) 
qs 


Introducing in Eq. (41) last two parts inside of exponential function, Eq. (26), we find the final 
expression of the function Ze, 


= [22h _ fe _ Sis Ph Py Ae 


Gq qt 
_ Vig + Yin — 2(Yej + ss) 
q3 


Our goal is to calculate the integral: 


Le ep ='({] {| AX, AY AX, | AY {AX AV KA 


—VriXkj + XK Kj 


oy a A OS 
(xe; + vie) ce + Yi i 


Using a formula below: 
+00 
| XN e~ PX? +24Xx dy = 
(3) 


0 


we find 

ba 3x2, 2( xn; + Xp1)xX 

de exp | — 3xk0 4 2%ks + Xi) *e0] _ 
ko€Xp 2 2 
—0o qi q1 
2 

[ex (xij + Xx) 

and 


ae 3yK 2(Vxj a Yur) Veo 
[ AV Ko VRoeXP - —2 as “ee pe i 


[7 
3 


2 
(Vej + Yet) (Vx; a Yet) 
Taking into account the expressions Eq. (45) and Eq. (46), the integral | a 2 transforms into form: 


3q3 3 
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TT 
Le ep = 3 1192 {|| | AX, AYR AX KI AV KL x 


—VeiXKj + XKVKj (Vx; ie Yet) 
(xz, + Ve; Oks + Yiu) 3 


with the new expression for Zeyp : 


exp 


2 
Zexp = XP soa hg hig) = 
1 
(yz; + Via — Veg Ve) | - 
3q3 "7 : 
The first term in the sum —YxjX_j + X41 Y_j Of the last expression for i does not depend on 


variable x;;. Therefore, we can calculate the integral 


— 1 204(,.2 

(XE - XR IX 

Leas = | Oi —— 3 ( kl—*kj xl) 
—00 Xi + Viet 

For this purpose, we use the definition of Gamma function ['(x) : 


+00 
mt) = | der tee 
a* 0 


to rewrite the [,,, in the form: 


1 +00 5 +00 
= = T 
Lis = ra |. dte kl ie AX; x 
204 204 
exp |- (= + t) xg, + "ae xe 


Using again a formula, Eq. (44), one obtains the result for I, : 


l= ur. (ie fy Se x 
TA) Jo J2a,/3+T 


exp (as) /(2a,/3+ D| 


3] 


In this expression for /,,,, at variation of variable T in the limits from 0 up to +00, the function 


exp (“24 /(2a,/3 + D| 


: 2 Aa : De eee 
will change from e1%iej/3 up, to | . However, at T > +00 limit, the function e~”*!* will make zero a 


whole integrand of I,,,. Therefore, one can assume 


AsXpir2 
exp ( =) /(2a,/3 +0] w eMrkj/3 


and thus take the approximate expression for integral 


2 
= +00 —VEIT 
TU e “kl 

carb | dt 
0 


“Gy. (2a,/3 +t 


Deng r(1) 
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2 a 2 
We take into account the exponential function e“**//? from Eq. (52) in the factor e~?%/3*ki, where 
a, = 1/q?, of the expression Zexp, and together with obtained this factor the integral over dx,; of 
Ve ep with the first term in the sum —yxjXqj + XK1VKj Bives: 


+00 

Xkj 2142, 
ae Xj + Vij 

and therefore, we find that / fon = 0. 


We consider the second term of the sum —YyxjX_j + XxYxj in the integrand of vy a It is easy to 
show that 


+00 
ees = | - tt 9 aati), 
d(2a,x,;/3) hes Xt + Vict 

Eq. (52), 


one obtains d/d(2a1X%j/3)lex) ~ xyes, Therefore, using Eq. (53), we find again Gs = 0. 


if expression for /,,, is taken from Eq. (49). However, from the final expression for /,,,, 


We demonstrated / ae = 0 at calculating the average quantities for bey = Yxo, related to the last 


term in the Hamiltonian H,. One can show that the same average quantities, however, calculating 
now for im = Xxo, give also 17 = 0. To get this result we used the expression 
, exp 


1 
Te ep 3 11.72 {| i AX, ; AV {AX KAYE X 
—VeiXkj + XV j (xx + Xia) 4 
(xt tye JGmt ye) 3 = 
and 
+00 ; es 
| Vij ve pen 3k =0 
7. HARMONIC OSCILLATOR EQUATION FOR THE CENTRE-OF-MASS. 


Substituting in the Schrédinger equation, Eq. (13), results of average quantities for fey calculated in 
the above three sections, we find the equations of motion for the xg coordinate 
2 


Xo ; 71 Cy 
i> — Arqixo SG = i> — GAai 1 
Ay 


5 xoq? — iv(N —1)C, + VIN - 1)C, 


and Yo coordinate 


q3 


Jy ' , {fC 
a A2q2¥0 — > & = — C,A2q3 + 
0) 2 7 
2 Yod2 — iW(N — 1)C, —V(N — 1)C, 
of the centre-of-mass. In equations Eqs. (57) - (58) A, = 4 and Az = 1/A 


Equating imaginary parts of both these equations, we find: 
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X 
i oe 1)C, 
Yo C2 
Es = > UN = 1)C,, 
from where 
Xo 
C = — 
b 
Yo 
C =— 
ss b 


Here, we introduced the constant b = 1 — 2v(N — 1). From Eq. (60), we express Xp and yp through 
the C, and C,, respectively, and substitute them in the real parts of Eq. (57) and Eq. (58). We obtain 


C 
C, + Aqxo = 2v(N — 1) (20.4, = =) 


qi 
’ Cy 
Cz + A2Vo = 2V(N — 1) (20.4, + 4). 
2 
Our goal is to consider solution of these equations on the first order small quantities, therefore, we 
omit the C,A, and CA, terms from the consideration and assume that q? = q?) and q% = qo. 
Taking into account the relationship, Eq. (60), we write a set of equations: 


x 2v(N — 1)y 
XO Fie? ao ( : )¥o 
b Tob 
Yo 2v(N — 1)Xo 
pee = ae 
9120 
We try to find the solutions in the form x9 ~ e'@¢ and yo ~ e'@* then Eqs. (62) reduce to 
—@* +4, = -iK,o 
—@* +2, =iK,6, 


where @ = w/b, A, = A,/b, Ay = Az/b, K, = 2v(N — 1) /(bq?y) and K = 2v(N — 1)/(bq3p). 
Multiplying two equations of Eq. (63) to each other, one obtains 
(—&? + 4,)(—6? + A,) = Ky, K26? 
and thus the equation 
@o* — (A, + a; + K,K) + Aaa = 0. 
The solution of this equation is: 
Less & 

(@*)12 =. 5 At Age K,K2) + 

sere & Oo yo 

F (A, + Az + K,K) i,’,| . 


We analyse an effect of the different cases of statistics of particles v and an harmonic potential 
anysotropy A, and A, on centre-of-mass oscillatory frequency w?. Let assume that we consider the 
system of bosons v = (. For this case of particle statistics, b = 1, K, = Kz = 0 and from equation 


(0*)12 > Wi 2 = 1/2(A, + Az) £1/2(A, — Ad), 
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i, =A; 
Oe 2s: 


For the case v = 0 and isotropic harmonic potential A, = A, = 1, we have w? = w3 = 1. 


For the system of anyons v # 0 and arbitrary harmonic potential 2, and A, the centre-of-mass 
oscillatory frequences are determined by Eq. (64). 


8. CONCLUSION 


So, 


after scrutinizing problems and tasks this work and with the help of acquired results the 


following statements can be done to conclude the work: 


> 


= 


OO play Do 


10. 


11. 
{2 
13. 


In order to perform the tasks set up in this research a new method for calculating problems has 
been utilized. It is a cumulant method to get the width equation and equation of the motion for 
the centre of mass to get the collective frequencies. This method gave an opportunity to avoid for 
writing an action and take a variation over it and also solve a differential equation of the second 
order which is an overwhelming task. Instead, we solved integrals of Gaussian type which is 
much easier than a variational approach. 


By utilizing the cumulant method, the expression for the centre-of-mass oscillatory frequencies 
for the system of anyons and arbitrary harmonic potential has been derived in the last chapter. 
Furthemore, special cases for the determined frequency have been considered with the 
parameters v and b, Kj, K>. 
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